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Abstract. Measurements of the diffractive structure function, of the proton at HERA are used 
to extract the partonic structure of the Pomeron. Regge Factorization is tested and is found to 
describe well the existing data within the selected kinematic range. The analysis is based on the 
next to leading order QCD evolution equations. The results obtained from various data sets are 
compared. An analysis of the uncertainties in determining the parton distributions is provided. The 
probability of diffraction is calculated using the obtained results. 
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INTRODUCTION 

In the last lOvears a large amount of diffractive data was accumulated at the HERA 
collider inlGlBn. There are three methods used at HERA to select diffractive events. One 
uses the Leading Proton Spectrometer (LPS) [jsj] to detect the scattered proton and by 
choosing the kinematic region where the scattered proton looses very little of its initial 
longitudinal energy, it ensures that the event was diffractive. A second method [ 2] simply 
requests a large rapidity gap (LRG) in the event and fits the data to contributions coming 
from Pomeron and Reggeon exchange. The third method ^ relies on the distribution 
of the mass of the hadronic system seen in the detector, Mx, to isolate diffractive events 
and makes use of the Forward Plug Calorimeter (FPC) to maximize the phase space 
coverage. We will refer to these three as ZEUS LPS, HI and ZUES FPC methods. 
The experiments 0, IH 0] provide sets of results for inclusive diffractive structure 

function, XpF2^^\ in different regions of phase space. In extracting the initial Pomeron 
parton distribution functions (pdfs), the data are fitted assuming the validity of Regge 
factorization. 

In the present study, Regge factorization is tested. New fits, based on a NLO QCD 
analysis, are provided and include the contribution of the longitudinal structure func- 
tion. The obtained PDFs are systematically analyzed. A comparison of the different 
experimental data sets is provided. Additional quantities derived from the fit results are 
also presented. 

In order to make sure that diffractive processes are selected, a cut of Xp < 0.01 was 
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performed, where Xp is the fraction of the proton momentum carried by the Pomeron. It 
was shown [7] that this cuts ensures the dominance of Pomeron exchange. In addition, a 
cut of 2^ > 3 GeV2 was performed on the exchanged photon virtuality for applying the 
NLO analysis. Finally, a cut on Mx > 2 GeV was used so as to exclude the light vector 
meson production. 



REGGE FACTORIZATION 



The Regge Factorization assumption can be reduced to the following. 



(1) 



where f{xp, t) represents the Pomeron flux which is assumed to be independent of /3 and 
and F{fi,Q^) represents the Pomeron structure and is /3 and dependent. In order 
to test this assumption, we check whether the flux f{xp,t) is indeed independent of /3 
and on the basis of the available experimental data. 

The flux is assumed to have a form ~ x^^ (after integrating over t which is not 
measured in the data) . A fit of this form to the data was performed in different 
intervals, for the whole /3 range, and for different /3 intervals for the whole range. 

Figure [l] shows the dependence of the exponent A for all three data sets, with the 
Xp and Mx cuts as described in the introduction. The HI and the LPS data show no 
dependence. The ZEUS FPC data show a small increase in A at the higher region. It 
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FIGURE 1. A as a function of g for Xp < 0.001 and Mx > 2 GeV, for the three data sets, as indicated 
in the figure. The hne corresponds to a fit over the whole region 

should be noted that while for the HI and LPS data, releasing the Xp cut to 0.03 seems 
to have no effect, the deviation of the ZEUS FPC data from a flat dependence increases 
from a 2.4 standard deviation (s.d.) to a 4.2 s.d. effect (not shown). 

The /3 dependence of A is shown in figure |2l All three data sets seem to show no /3 
dependence, within the errors of the data. Note however, that by releasing the Xp cut to 
higher values, a strong dependence of the flux on /3 is observed (not shown). 

We thus conclude that for jc^ < 0.01, the Pomeron flux seems to be independent of 
and of /3 and thus the Regge factorization hypothesis holds. 



NLO QCD FITS 

We parameterized the parton distribution functions of the Pomeron at = 3 GeV^ in a 
simple form of Ajc*(1 —xY for u and d quarks (and anti-quarks) and set all other quarks 
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FIGURE 2. A as a function of /3 for Xp < 0.001 and Mx > 2 GeV, for the three data sets, as indicated 
in the figure. The fine corresponds to a fit over the whole j3 region 



to zero at the initial scale. The gluon distribution was also assumed to have the same 
mathematical form. We thus had 3 parameters for quarks, 3 for gluons and an additional 
parameter for the flux, expressed in terms of the Pomeron intercept ap(0). Each data set 
was fitted to 7 parameters and a good fit was achieved for each. The HI and ZEUS FPC 
had X^/^i ~ 1 J while for the LPS data, the obtained value was 0.5. The data together with 
the results of the fits are shown in figure|3l The following values were obtained for ap(0). 




FIGURE 3. The diffractive reduced cross section of the proton multiplies by Xp, as a function of Xp for 
the different data sets (the most right plot is for the LPS data) in different bins of and j3, as indicated 
in the figure. The bands are the results of the fits including uncertainties. 

for each of the three data sets: ap(0) = 1 . 138 ± 0.01 1 , for the ZEUS FPC data, ap{0) = 
1.1 89 ±0.020, for the ZEUS LPS data,aj,(0) = 1.178 ±0.007, for the HI data. 

The parton distribution functions are shown in figure IH for the HI and the ZEUS FPC 
data points. Because of the limited /3 range covered by the LPS data, the resulting pdfs 
uncertainties are large and are not shown here. For the HI fit one sees the dominance 
of the gluons in all the /3 range. For the ZEUS FPC data, the quark constituent of the 
Pomeron dominates at high /3 while gluons dominate at low /3. We can quantify this by 
calculating the Pomeron momentum carried by the gluons. Using the fit results one gets 
for the HI data 80-90%. while for the ZEUS FPC data, 55-65%. 



FIGURE 4. Quark and gluon pdfs of the Pomeron as obtained from the HI data fit (left two figures) 
and from the ZEUS FPC data fit (two rightmost figures) as a function of j3, at different values of . 

COMPARISON OF THE DATA SETS 

One way of checking the compatibility of all three data sets is to make an overall fit for 
the whole data sample. Since the coverage of the /3 range in the LPS data is limited, we 
compare only the HI and the ZEUS PC data. A fit with a relative overall scaling factor 
of the two data sets failed. Using the fit results of one data sets superimposed on the 
other shows that the fit can describe some kinematic regions, while failing in other bins. 
This leads to the conclusion that there seems to be some incompatibility between the 
two data sets. 

PROBABILITY OF DIFFRACTION 

It is of interest to calculate the probability that a certain parton is produced in a diffractive 
process [8]. The probability of diffraction on quarks and gluons, as a function of Bjorken 
X at different values of 2^ are shown in figure|5l using the results of the HI and the ZEUS 
FPC data fits. The ZEUS FPC data shows that throughout the whole kinematic range 
shown in the figures, the probability for diffraction is not bigger than 0.15, far from the 
Pumplin [9] limit of 0.5. This is not the case for the HI data for which, at small x and 
low 2^, the probability of diffraction induced by gluons becomes greater than 0.5 and 
thus unphysical. Note however that the results for x <2- 10^"^ are in a region where HI 
has no data and thus the calculated probability in this region is an extrapolation based 
on the fit parameters. In order to get physical results, some process, like saturation, must 
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FIGURE 5. Probability of diffraction as a function of x, at different values of Q^, calculated from the 
results of the HI data fit (left figure) and from those of the ZEUS FPC data fit (right figure). 



lower the expected value. 
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